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ABSTRACT 

After two ALMA observing cycles, only a handful of [C ii] 158 /im emission line searches in z > 6 
galaxies have reported a positive detection, questioning the applicability of t he loc al [C ii]-SFR relation 
to high-z systems. To investigate this issue we use the Vallini et al.| (120131 V13) model, based 
on high-resolution, radiative transfer cosmological simulations to predict the [C ii] emission from 
the interstellar medium of a z « 7 (halo mass = 1.17 x 10^^ Mq) galaxy. We improve the 
V13 model by including (a) a physically-motivated metallicity {Z) distribution of the gas, (b) the 
contribution of Photo-Dissociation Regions (PDRs), (c) the effects of Cosmic Microwave Background 
on the [C ii] line luminosity. We study the relative contribution of diffuse neutral gas to the total 
[C ii] emission (Ediff/Ftot) for different SFR and Z values. We find that the [C ii] emission arises 
predominantly from PDRs: regardless of the galaxy properties, Fdig/Ftot < 10% since, at these early 
epochs, the CMB temperature approaches the spin temperature of the [C ii] transition in the cold 
neutral medium (Tomb ~ ^ 20 K). Our model predicts a high-z [C ii]-SFR relation consistent 

with observations of local dwarf galaxies (0.02 < Z/Zq < 0.5). The [C ii] deficit suggested by actual 
data (Lcii < 2.0 x lO^L© in BDF3299 at z « 7.1) if confirmed by deeper ALMA observations, 
can be ascribed to negative stellar feedback disrupting molecular clouds around star formation sites. 
The deviation from the local [C ii]-SFR would then imply a modified Kennicutt-Schmidt relation in 
z > 6 galaxies. Alternatively/in addition, the deficit might be explained by low gas metallicities 
(Z < 0.1 Zq). 

Subject headings: galaxiesihigh-redshift, galaxiesiism, cosmology:theory, submillimeteriism, 
line:formation, cosmology:observations 


1. INTRODUCTION 

The study and characterization of the interstellar 
medium (ISM) of galaxies that formed in the early Uni¬ 
verse is entering a golden era thanks to the unprece¬ 
dented capabilities of the Atacama Large Millimeter- 
submillimeter Array (ALMA). In particular, the 158 /im 
emission line due to the ^^3/2 — P 1/2 fine-structure 
transition of ionized carbon ([G ii] ), being the domin ant 
coolant of the neutral diffuse ISM dWolfire et al .|200^, is 

a (IStac 


by far the b rightest line in the far-intrared band (Stacey 
et al. 1991). In addition to the diffuse neutral gas, the 
[G iij line can be excited in other components of the in¬ 
terstellar medium such as high density photodissociation 
regions (PDRs), and in the diffuse ionized gas, where the 
main driver of the [C ii] emissivity are the collisions with 
free e~. Although precisely assessing the relative contri¬ 
bution of the various gas phases to the total line emission 
might be difficult, [C ii] line remains an exquisite proxy 
to characterize the ISM of galaxies that form ed during 
the Epoch of Reionization (EoR; z « 6 —7) (e.g Carilli & 
Walter 2013). Before the ALMA advent, the [G iij line 
from z > 4 was solely detected in galaxies with ext reme 
star formation rates (SFRs) («1000 M0yr“^) (e.g. Gox 


^ Scunla, Nnrma.le Siineriore. Pisa,. Tla,lv 

^ Kavii IFM U (W FI), Todai Institutes tor Advanced Study! 

tlie University ot ibkyo 


et ai.||2011| 

Carilli et al.||20131 |Carniani et al.||20131 |De| 

Breuck et a 

. 2014), or in those hosting 

Active Galactic 

Nuclei (AGN) (e.g. Maiolino et al. 12005 

IVenemans et al. 

2012 Gallerani et al. 2012| Gicone et a 

|2015). 

In the tirst years ol ALMA operations, the [CJ iij emis- 


sion line has been detected in a handful of galaxies 
with modest star formation rates (50 — 300 Mq yr“^) at 


EoR 

Garilli et al. 

2013 Garniani et al. 

2013 

et al. 

2014 

Riechers et al.|2014). Viceversa, 


forming galaxies (NSFGs; SFR wlO Mf7,yr~^) at the end 
of Epoch of Reionization (z > 6) (e.g. Walter et al.|2012 
Knnelgj_gt ^y|20j_3RGpnzaiezHjC)pez_eUai.J^^ |Uuchi| 
et al 112013 Ota et al.||2014 Schaerer et al. 2015). These 


early results seemed to be at odds with the correlation 
between the intensity o f the [C ii] line and the SF R ob¬ 
served in local galaxies (De Looze et al. ]MI||M4l ), thus 
questioning its applicability to sources at z ^ 6. Only 
very recently, three different ALMA campaigns target¬ 
ing z « 5 — 7 Lyman Alpha Emitters (LAEs) and Ly¬ 
man Break Galaxies (LBG s) have yielded [C ii] detec¬ 
tions: Maiolino et al. ( 2015|) in th e vicinity of BDF3299 


a LAE at z « 7.1, Gapak et al. J201^|in a sample of 


LAEs at 5.1 < z < 5.7, and Willott et al. (2015) in two 
luminous LBGs at z « 6, being m agreernent with the 
[C ii] luminosity expected from lower-z observations in 
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star forming galaxies. 

In the nearby Universe, the [C ii]-SFR relation holds 
for a wide range of galaxy types, ranging from metal 


frared galaxies, and AGN hosting galaxies (Boselli et al. 

2002[ [De Looze et al.[ [20111 jSargsyan et al.| [20121 |Ue 

Looze et al.| 20141 

Pineda et al.l [20141 

Herrera-Camus 

et al.|[2015). TEeTC 

ii] emission from PURs is primarily 


OB stars tha t form in the vicinity of the photodissocia¬ 
tion regions (|Hollenbach & Tielens 1999). The relation 
between SFR and the [O iij luminosity in the neutral 
diffuse gas is more subtle. On the one hand, the [C ii] 
emissivity is proportional to the gas heating due to the 


sitv of the EUV radiation ( 

Wolfire et al. 2003 

Herrera- 

Camus et al. 

2L)15I. Un the other hand, an increasing 


neutral medium (CNM) with respect to the warm neu- 
tral medium (WNM ), thus reducing the [C ii] luminosity 
Observational studies have found 
of the Galaxy, t he [C ii] emission i s 


[2M^. 

ilaneoi 


(Vallini et al. 

that, in the p „ , _ 

mostly associated with dense PDRs (|Pineda et al.|2013|) . 
On the contr ary, in low metallicit y local dwart galaxies 


(e.g. Haro 11, Cormier et al.|20l'2 |, nearby galaxies (e.g. 
M51 and M3 1, IKramer et al. | |2(J13[ [Parkin et al.|[20T3| 


Kapala et al. 


26l5|and the outskirts of the Milky Way 


(jLanger et al 2(J14[) the PDR contribution can be as 


small as « ll)%. 'I'he [C ii]-SFR relation in these cases is 
shallo wer than that of starburst galaxies (De Looze et al. 

Mil. 


From a theoretical point of view the intensity of the 
[C ii] line from high-z galaxies ha s been computed both 
through numerical simulations ( e.gjNagamine et al.|2006|) 
and semi-analytica l models (e.g. Gong e t al.|2012| I Vluhoz 
&: Furlanett o|2014l [Popping et al.|2U14|). Recently loisen 
et al.| ([2U15|) present a multi-phased ISM model consist- 
ing 01 molecular clouds, embedded within a cold neutral 
medium of atomic gas, and hot, partly ionized gas. The 
model, applied on top of a cosmological SPH simulation 
of massive star-forming galaxies on the main-sequence at 
z = 2, self-consistently calculates the relative contribu¬ 
tion of the various phases to the [C ii] emission. 


In the previous paper of this series Vallini et al. (2013) 
(hereafter, V13) computed the [C iij emission arising 
from the neutral diffuse gas of a single prototypical high- 
z {z ~ 6 . 6 ) galaxy, extracted from a SPH cosmological 
simulation, further implemented with radiative trans¬ 
fer calculation. This is crucial to model the intensity 
of the galaxy internal UV field and the consequent gas 
ionization structure. The calculation of the [G ii] emis¬ 
sion is performed thanks to a sub-grid model describing 
the thermal equilibrium of the cold and warm neutral 
medium as a function of the FUV radiation field inten¬ 
sity within the galaxy. The spatial resolution (« 60 pc) 
allows to properly describe the ISM small-scale density 
structure. Here we present an updated version of the V13 
model that allows us to also compute the [C ii] emission 
arising from the clumpy molecular gas, and the effect of 
the increased GMB temperature on the [C ii] observabil¬ 
ity. The aim is to finally assess whether the local [G iij- 
SFR relation holds at high-z, and what we can learn from 
any deviation from it. 


2. MODELLING [C ii] EMISSION 

In this Section, we first summarize the main character¬ 
istics of the V13 model, referring the interested reader to 
V13 for further details. Next, we describe the additional 
features implemented by this work. 


2 . 1 . The VIS model 


We run GADGET-2 (Springel 20051 cosmological SPH 
hydrodynamic simulations of a (10/i“^Mpc)^ comoving 
volume with a mass resolution of 1.32 ( 6 . 68 )xlO^M 0 
for baryons (dark matter). We take a snapshot at red- 
shift z = 6 . 6 , identify the most massive halo (total mass 
Mh = 1.17 X 10^^ Mq, Tyir « 20 kpc), and we select 
a (0.625 /i“^Mpc)^ comoving volume around the center 
of the halo. We post-process our simulations with the 
UV radiative transfer (RT) code LIGORIGE. Gas prop¬ 
erties are resolved on a fixed grid with a resolution of 
^ 60 pc. We complement the simulation with a sub-grid 
model taking into account the cooling and heating pro- 
cesses producing a multi-p hase thermal ISM structure 
([Wolfire et al.[ [19951 |2003[) . According to this model, 
the neutral gas in the ISM is constituted by a two-phase 
medium in which the cold neutral medium (CNM) and 
the warm neutral medium (WNM) are in pressure equi¬ 
librium. The relative abundance of these two compo¬ 
nents depends on (i) the gas metallicity, Z, determining 
the coolants abundance, and (ii) the FUV flux, Gq, in the 
Habing (6 — 13.6 eV) band, controlling the photoelectric 
heating produced by dust grains. The value of Gq scales 
with SFR and is calculated as: 


rl3.6eV , 

Go(r, SFR) = SFR x ' 

47 r|r —r 


,:di/ 


( 1 ) 


where n* is the number of sources, fl is the positions, l,^^i 
is the monochromatic luminosity per source. We com- 
pute l ,^_i by using STARBURST99 template Leitherer et al. 
(1999), assuming continuous star formation fST’), an age 
lOMyr for the stellar populatior 0 and setti ng th e metal¬ 
licity accordingly to the cell value (see Sec 2.2). In this 
work we explore the range SFR= [0.1 — 100] Mq yr“^. 
For each value of the star formation rate and metallicity 
we run the V13 sub-grid model computing the expected 
distribution of the cold and warm diffuse gas within the 
galaxy. 

We calculate the [C ii] emissivity through Eq. (3) of 
V13. We note that this equation is valid under the as¬ 
sumption that the density of the colliding species (elec¬ 
trons and hydrogen atoms) is much lower than their cor¬ 
responding critical densitie^ and that no external ra¬ 
diation field is present. While the first assumption is 
justified by the fact that in our simulation rie <C = 


^ We keep fixed the original assumption of V13 in which a con- 
tinuous SF with an age of £* = 10 Myr for the stellar p opula¬ 
tion was ba s ed on the Himiko SED fitting presented in |Ouchi| 
et al. I (12009112013 ||. T his scenario have been recently connrmed 
Zabl et ai.| and |Schaerer et al.| (|2015|| that, when adopt- 

ing continuous/exponential rising (decliningj star formation histo¬ 
ries for Himiko, obtain age of th e stellar population in the range 
= 1 — 35 Myr l |Zabl et al.|201^ , and £ * = 10 — 40 Myr l |Schaer^ 

‘i^he critical density for collision with neutral hydrogen atoms 
(or with e~) has been computed at T = 100 K, a value consistent 
with the CNM temperature. 
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Fig. 1.— Panel (a)-. The black solid line and the gray shaded 
region represent the radial profile of the baryonic overdensity (A) 
and its r.m.s. fluctuation for our simulated galaxy. The distance 
is rescaled with the virial radius of the galaxy. With lig ht cyan 
lines we p l ot the m ean profiles for the ~ fOO galaxies in t helPallot- 
tini et ahj l| 2014a| (P14)) simulation (see also Fig. 1 in dPallottim 
et al.||i2Ul4bT rr' "we highlight with solid dashed/dot-dastied lines 
two most massive P14 galaxies. Such galaxy are hosted in dark 
matter halos with mass log(M;i/M 0 ) 12 and log(M^/M 0 ) 11 
respectively, comparable the one of the galaxy adopted in this work 
(log(M/j,/M 0 ) PS 11.1). Panel (b): Metallicity {Z) radial profile for 
the P005’ model. For the P005, Z is calculated using the density 
of our simulated galaxy and rescaling Z — A relation found in P14. 
The relation is rescaled so the mean metallicity of the galaxy is 
{Z) = 0.05 Z 0 (grey dashed line). Note that (Z) is calculated con¬ 
sidering only those cells whose baryonic overdensity is A > 200, 
i.e. up to rjryir ^ 0.5. The MCs are, on average, clusterized near 
the center, and this affects their mean metallicity. 

8cm“^ and nn ^ 3000cm“^, in Sec. 2.3 we investigate 
whether the CMB may affect the intensity of the [C ii] 
emission. As pointed out in V13, the CNM accounts for 
« 95% of the total [C ii] emission arising from the diffuse 
neutral medium. Given this result, we refer to the diffuse 
medium as CNM. 


2.2. Metallicity 

The V13 model assumes that metals are uniformly dis¬ 
tributed within the galaxy. Hereafter, we refer to models 
with a uniform metallicity distribution as ^ “C-models”. 
In this work, we also consider the possibility that the 
distribution of metals follows the density distribution, 
by relating Z to the baryonic overdensity of each cell 
A = pg 

as I Pc{z), where pgas is the gas density in the cell 
and Pc{z) is the critical density at redshift z. This is in 
agreement w ith the parametrization adopted by [Keating] 
et al. (2014) to describe the circumgalactic mediu m of 
high-z galax i es, and with the results presented in jPal-j 
lottini et al. (2014a, P14) that has been applied to our 
galaxy. 

P14 use a customized version of the adaptive mesh re- 


finement code RAMSES ( 

Teyssier 

2002) 

vestigate the metal enric 

ament 0 

high- 


P14 star formation is included via subgrid prescriptions, 
and supernovae feedback is accounted by implementing 
a metal-dependent parameterization of stellar yields and 
return fractions based on population synthesized models. 


® For example, we indicate with C02 a model in which a uniform 
Z = 0.2 Z 0 is imposed. 


The P14 galaxy sam ple reproduces the ob s erved cosmic 
star f ormation rate ([Bouwens et al. 2012[ Zheng et al. 


2012) and stellar mass densities ( [Gonzaiez et al. 2011[ ) 
evolution in the redshift range A < z ^ 10. In the upper 
panel of Fig. [T] we plot the radial profile of the bary¬ 
onic overdensity (A) and its r.m.s. fluctuation of our 
simulated galaxy. We test the profile with a sample of 
^ 100 galax ies extracted from P14 (see also Pallottinij 
et al. 2014b[ in particular Fig. 1, upper panel). Among 
the P14 galaxies, with solid dashed/dot-dashed lines we 
highlight two galaxy hosted in dark matter halo with 
mass \og{Mh/MQ) « 12 and \og{Mh/MQ) « 11, respec¬ 
tively. This is comparable to the dark matter halo mass 
of the galaxy adopted in this work {log{Mh/MQ) « 11.1). 

In particular, P14 found a tight correlation between Z 
and A for log A .li 2, namely for overdensities typical of 
galaxy outskirts/ISM. While in the IGM (log A < 2.3) 
the metallicity is only weakly correlated with A, in the 
ISM (2.3 < log A < 4.5) the Z — A relation is tight. This 
is due to the fact that the most overdense regions denote 
the location in which stars form, and that are therefore 
more efficiently polluted with metals. We fit the Z — A 
relation provided in their paper and we normalize the 
relation to the mean metallicity (Z) over the galaxy, i.e. 
over those cells whose overdensity is A > 200. Hereafter, 
we refer to models that take into account this Z — A re¬ 
lation as “P-models”. The density-dependent metallicity 
case with {Z) = 0.05^0 is called P005, and the prohle is 
shown in the (b) panel of Fig. 

2.3. Molecular clouds and PDRs 

Beside the emission arising from the diffuse neutral 


todissociation regions 

Tielens & Hollenbach 

1985 

Hoi- 

lenbach & Tielens|1999 

around molecular clouds (i\ 


lb establish whether t 

he gas in a cell becomes gravita- 


tionally bound, we apply the Jeans instability criterium. 
We define molecular cells those satisfying the following 
condition: Mceii > Mj {Tcmb, riceii) , where Mceii is the 
mass of the gas in the cell, Mj is the Jeans mass at 
the CMB temperature (Tcmb) at z = 6.6, and at the 
density of the gas in the simulation cell (jiceii)- As a 
caveat we note that is possible that other processes such 
as the cosmic rays heating (e.g. Papadopoulos||2010 ) can 


increase the temperature above the CMB floor, even at 
z « 6. By applying this prescription, we find that the to¬ 
tal mass of molecular hydrogen in the simulated galaxy is 
Mh 2 = 3.9 X 1O®M0, consistently wit h previous theoret- 


ical estimates (Mh 2 = 2.5 x 10® M 0 , Vallini et al.||2012 | 
and observationa l constraints (Mh 2 < 4.9 x 10'^ M 0 


Wagg et a Lj[2 009[) on z « 7 LAEs. 
in 


Since we hnd that on average the simulation cells con¬ 
tain a molecular hydrogen mass ^ 3 x 10® M 0 , 

we consider each molecular cell as a Giant Molecular 


Cloud (CMC; A/gmc = 10® - 10® M 0 , Murray 2011). 
The properties of GMCs are controlled by a turbulerit 


mal shock waves (Padoan 

1995 

Ostriker et al. 

2001 

Padoan et al. 

2014 

). The problem of turbulent frag- 


cally 

(e.g. Krumholz & McKee|2005 

Padoan & Nordlund 

2011 

IHennebelle & Chabrier 

2mT' 

2U13D or numericaliy 

(e-g- 

Vazquez-Semadeni|1994 

Kim & Ustriker|2002 

Kim 
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Fig. 2.— Slices of thickness 57 pc cut through the simulated 
galaxy for the P005 model showing the metallicity distribution 
(color scale) and molecular clouds (black regions). Note that the 
metallicity of the overdense regions in which MCs reside is greater 
than the mean value (Z) = 0.05 Zq. 


et al. 2003 Wada 2008 Tasker & Tan 2009). Analyti- 
cal models as well as numerical simulations show that the 
distribution of the gas density (nd) in an isothermal, non 
self-gravitati ng, turbulent medium follows a log-norm al 
distribution ( Padoan||1995 |Padoan fc NordlundpOlI ): 


p(lna;)c?lna; = 


(27rcr2)i/2 


exp 


In a: — In a: 


_ ( 2 ) 

where x = rid/rid, rid — ri^M/^, uq is the average 
number density of the CNM (uq = 50cm“^, se e V13), 
A4 = 10 is the Mach number value suggested by Kainu- 
lainen & Tan (2013), the turbulent velocity dispersion 
is given by cr = A/ln[l -|- with /3 = 0.5 and 

Ina; = —O.Str^. If we assume that each GMC in our sim¬ 
ulation is composed by a set of clumps, we can compute 
the densities rid of each clump by adopting an iterative 
approach that consists of the following steps: 


1. Select rid by sampling the density distribution (Eq. 

2. Set the clump radius equal to the Jeans length rd = 
^j{Tcmb, rid)- 


3. Calculate the clump mass Md = {9 ,tt/‘ i)mHndr/.i. 

4. If Md < Mceii calculate the residual mass in the 

cell = Mceii - Md] go to step 1 . 

5. If Md > Mceii reject the density sampled and as¬ 
sume Md = Mceii- Calculate rd as in step 2. and 
use it to derive the resulting clump density. 


Through this procedure, we find that the MCs in the 
simulations are characterized by the following properties: 
(logKi/cTO-3)) = 2.9 ±0.4, (Md) = 50±20Mo and 
radius {rci) = 0.7 ± 0.3 pc. 

In Fig. 0 we show four 57 pc-thick slices of the simu¬ 
lated box lowing the metallicity distribution and molec¬ 
ular clouds (black regions) for the P005 model. By in¬ 
specting the maps shown in the Figure we note that 


molecular cells reside in highly overdense regions mainly 
clustered at the center of the galaxy. The molecular cells 
are located predominantly in the innermost region of the 
galaxy {d-MCs ~ Ikpc) and have Z k- 0.2 Zq, whereas 
the mean metallicity over the galaxy is {Z) — 0.05 Zq. 
This implies that any line arising from the PDRs would 
provide an upper limit on the mean metallicity of the 
galaxy. Finally, to calculate the [C ii] emis sion from 
PDRs, we couple our simulati on with UCL_PDR (Bell et al. 
2005 2007 Bayet et al.|2009) a PDR code that allows to 
derive the [C iij emissivity as a function of the intensity 
of the FUV radiation field (Gq, Sec. 2.1), metallicity {Z, 
Sec. 2.2) and molecular gas density {rid, Sec 2.3). 

2.4. CMB effects on [C ii] emission 

Since the CMB sets the minimum temperature of the 
ISM to Tcmb = 7cmb(1 + s-f high z it represents 
a strong bac kground against whic h the line fluxes are 
detected (e.g. da Cunha et al.|2(iT3 ). The contrast of the 
cloud emission against the CMB radiation is given by the 
following relation: 

A4 = [R,(r,) - B,(rcMB)] (1 - e-"-'), (3) 

where Tg is the excitation (or spin) temperature. If we 
assume that the [C ii] line is optically thin in the sub¬ 
millimeter, i.e. « 1 — Tj,, the ratio ((1 ± z)Fb ~ 

Iv/d\)) between the flux observed against the CMB and 
the intrinsic flux emitted by the cloud is (see also |da| 
Cunha et al.||2013 ): 

-Bi/(7cmb) 


[B.jTe) - B,{Tcmb)]t, 
^ B,{T,)t, 


= 1 - 


B,{Ts) 


ri-i B-c^iIiy ± ± rieC , ± 


''ul 


Eq. (4|) clearly shows that as the Tcmb approaches Tg 
the observed flux tends to zero. The ratio of the pop¬ 
ulation of the upper (^± 3 / 2 , labeled with u) and lower 
{’^Pi/ 2 , labeled with 1) levels of the [C ii] 158 /im transi¬ 
tion formally defines the spin temperature: 

riu _ Biulv ± UeCf^ ± rinC// _ Qu 

^ u — ^ 5 

9 i 

where A^i is the Einstein coefficient for spontaneous 
emission, B^i {Biu) is the stimulated emission (absorp¬ 
tion) coefficient, Gf^ (G^) is the collisional excitation 
rate for collision with e~ (protons), G®; (G^) is the 
collisional de-excitation rate for collision with e~ (pro¬ 
tons), and Ue {nn) is the number density of e~ (pro¬ 
tons). For the [C ii] 158 ^m line emissionA^; = 2.36 x 
10-6s-i, GfjT) = (8.63 X 10-6/ff;Vr)7iJr)e-^-/^, 


with 7 i„(T) « 1.6 if 1 00 < r < 10 ^ dCong et al.||2012| ) 
C//{T) is tabulated in Dalgarno & McCray (1972), and 
T* = huui/h « 91 K. 

In local thermal equilibrium (LTE) the collisional exci¬ 
tation and de-excitation rates are related by the following 
expression that depends on the kinetic temperature T : 


^e,H _ 9 u 

'^lu — „ ^ '-'ul ■ 

9 i 


By combining the equations ([^ and ([^ we obtain: 

A„i(l ± 1^) ± ^^eG®; ± UhC/ 


— = In — 
Ts A 


-Iff 

-'ul 


ui{$^)+n. 








C/^,e-¥ 


( 6 ) 


(7) 
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Fig. 3.— [C ll] spectrum for the P005 model assuming SFR = 
lM 0 yr“^ (upper panel) and SFR=lOM 0 yr“^ (lower panel) and 
rebinned over 5 kms“^ velocity channels. The emission from PDRs 
(diffuse neutral medium) is plotted in light red (dark blue). 


As discussed by Gong et al. (2012), the soft UV back¬ 


ground at 1330 A (/uv) produced by stars can in prin¬ 
ciple pump the [C ii] ions from the energy level ^f’ 1/2 
to This pumping effect, can lead to the [C ii] 

fine structure transition ^D ^/2 P 3/2 — Pi/ 2 ^ which 

would mix the levels of the [C ii] 158 /tm line and thus 
modify the Eq. ([^. However, the UV pumping effects 
are negligible in our calculations since the UV intensity 
inside the galaxy for all the SFR values considered is 
much smaller than the critical value for this effect to be- 
come important, namely 10 “^®erg cm“^Hz“^ sr“^ 
( Gong et al.||2012 ). 

We calculate the spin temperature of the [C ii] transi¬ 
tion in the PDRs by substituting into Eq. 0 
and T as resulting from the UCL_PDR outputs. The gas 
temperature within PDRs depends on the radius con¬ 
sidered and on the SFR and ranges between 20.7 < 
(T) < 800 K. We find that - 30 - 120 AT for 

SFR= 0.1 - lOOMayr-i. The in the GNM is cal¬ 
culated by considering the Ue, nn, and Tj. provided by 
the V13 sub-grid model. We obtain ^ 22 — 23 K 

approximately constant in the range of SFR considered. 
Since at z « 6.6 Tcmb ~ 20.7 K, we find that the [G ii] 
emission arising from PDRs is only slightly affected by 
the CMB (C « 0.8—1.0). Viceversa, the GNM is strongly 
attenuated at this redshift (C « 0.1 — 0 . 2 ); in this case 
the GMB effect becomes negligible only for galaxies at 
z < 4.5. 

3. RESULTS 

In Fig. 1^ we show the [G ii] spectrum obtained from 
the P005 model, assuming SFR = 1 Mq yr“^ (top panel) 
and SFR = lOM 0 yr“^ (bottom panel). In this Fig¬ 
ure, the contribution to the [G ii] emission arising from 
PDRs and the GNM is shown in light red and dark blue, 
respectively. The emission from PDRs arises predomi¬ 
nantly from the center of the galaxy, covering the ve¬ 
locity channels around ~ 0kms“^. The second peak 
in the PDR emission at u ~ 200kms“^ is produced by 
MGs located in the GNM clumps at the periphery of 
the galaxy (see Fig. [^. [G ii] emission from the diffuse 
medium, visible as the two peaks around ^ 100 kms“^ 
and ^ 200 kms“^ is instead always displaced from the 



Fig. 4.— Fraction of the [C ii] flux arising from the diffuse 
medium (Fdiff) over the total flux (Ftot), for four different models, 
as a function of the SFR, with (dark blue) and without (light blue) 
taking into account CMB effects on the [C ii] emission. 

center of the galaxy. The [G ii] line is relative^ narrow, 
with a FWHM~ 50kms“^, as in V13. In Fig. we plot 
the relative contribution of the diffuse medium to the to¬ 
tal [C ii] emission, Fdifi/Ftot, as a function of the SFR, 
for different G- and P-models, taking into account the 
GMB attenuation (dark blue). We find that the [C ii] 
emission in z « 6 galaxies is dominated by PDRs, since 
the GNM contribution is always < 10%, regardless of the 
metallicity profile and SFR considered. 

When the CMB attenuation of the GNM luminosity is 
negligible (i.e . typically for sources located at z < 4.5 
see Sec. 2.4), we find Fdis/Ftot = 0.05 — 0.45, con- 


nearby galaxies (ICormier et al. 120121 

Kramer et al. 12013 

Parkin et al.||2013[ Langer et al.||201^ 

5 . 


3.1. The [C u]-SFR relation 

In the previous Section, we have found that the [C ii] 
emission is dominated by PDRs, implying Lcii oc Mhj. 
We rescale the [C ii] luminosity of our “fiducial” model 
(SFR= lOM 0 yr“^, Mh 2 = 4 x 10®M 0 ) to an arbitrary 
molecular con tent by assuming the Kennicutt-Schmidt 
(KS) relation (Kennicutt et al.(l998 Kennicutt & Evans 


2012), namely a power-law correlation between the SFR 
and molecular gas surface densities, SsFR oc The 

range in power-law index (N) relating ^sfr and Eh 2 
depends on a variety of factors, among which the most 
important ones are the observe d scale, and the calibra¬ 


tion of star formation rates. Kennicutt et al. (1998) 


and Narayanan et al. (2012) report s uper-linear indeces 
N = 1.4 and fV « 2 respectively, while Bigiel et al. (2008) 


inferred an a pproximately l inear molecular KS relation. 
More recently Shetty et al. (2013), by performing a hier- 
ar chical Bayesian anal ysis on the same sample considered 
by |Bigiel et al.| (|2008|) conclude that N = 0.84, with 2cr 
raiige [U.63-1.UJ. In what follows we scale the molecular 
mass with the SFR by adopting N = 1 leaving to the 
last Section the discussion on the impact of different N 
on [G ii]-SFR relation. 

In Fig. we show the result of this procedure for dif¬ 
ferent metallicity profiles. Models with uniform metal- 
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Fig. 5.— [C ii] luminosities in solar units as a function of the SFR. Results from this work are shown with big diamonds and lines color 
coded as a function of the metallicity. Solid lines represent the result obtained assuming a constant metallicity: black for Z = Zq (Cl), 
orange for Z = 0.2 Zq (C02), and blue for Z = 0.05 Zq (C005). The magenta dashed line corresponds to a density-dependent metallicity 
(see Sec. 2.2) with (Z) = 0.05 Zq (P005). Model predictions are compared with data of local dwarf galaxies individually denoted by circles 
color-coded according to their metallicity (|De Looze et aL|2014||. The Icr scatter around the best fit relation for dwarf and local starburst 
galaxies is plotted in dark gray and light gray, respectively. Ureen empty (filled) points represent upper limits ('detec tions) of [C ill in LAE s 
and LBGs at 2 : ps 5 — 7. More precisely: the filled square represe nts the [C ii] detection nearby BDF3299 : 
filled hexagons show recent observations by 


limits on BDF3299, BDF512. and SDF46 (Maiolmo et a 


(|Ouchi et al.|2013 fota et al.i2 014|| (empty 
Ai7U3-zJJi and zo-Vii\ 11-5290 dbchaerer et 


Capak et al.|l^l5[l. filled triangles with red edges denote data 


H ar' 

al. 


repor ted bvIMaiolino et al .l l|2015|| , 
ita by |Willott et al!] | |2Ui5[ |. Upper 


are plot ted with empty squares, along with ttie uppe r limits on Himiko 


left triangle). 


star). iUK-i (lUta 
-ViiML)-529o ( jbcnaerer et al.|201^ (empty down-triang 


|2014[ | (up-triangle), HCM6A ( |Kaneka ret^^aJj 2013[l (empty di amond), 
iriangie), SDF1543 and SDF30ob i |Uonzalez-Li6pez et al.|2014[ | (empty 


licity are shown with a black dotted line for Z = Zq 
(Cl), orange solid line for Z = 0.2Zq (C02), and blue 
dotted-dashed line for Z = 0.05 Zq (C005). The results 
from our C-models are well described by the following 
best-fitting formula: 

logLcii = 7.0 -F 1.2 log(SFR) -F 0.021 log(Z) -F 

0.012 log(SFR) \og{Z) - 0.74 \og^{Z), (8) 


where Lqii is expressed in solar units, and the SFR in 
MQyr“^. The magenta dashed line indicates the pre¬ 
dictions for the Z-A relation with (Z) = 0.05Zq (P005). 
The slope of the [C ii]-SFR relation does not depend nei¬ 
ther on (Z) nor on the metallicity distribution. More¬ 
over, the [C ii] luminosity predicted by the POOS model 
is almost coincident with that obtained from the C02 
model, in the entire range of SFR considered. This 
can be understood by considering that, in the case of 
a overdensity-dependent metallicity profile the photodis¬ 
sociation regions located ad the center and dominating 
the emission, have Zpdr > (Z) and, more precisely, 
ZpDR « 0.2 Zq, namely the metallicity value of the C02 
model. 

We compare our predictions wi th the [C ii]-SFR rela- 
tions and their Itr scatter found by De Looze et al. (20141 
for local dwarf galaxies (dark gray shaded region) and lo- 
cal starburst g alaxies (light gray sha ded region). Dwarf 
galaxy data by De Looze et al. (2014) are shown through 
filled circles, ihdividually color-coded according to their 
Z. Upper limits on the [C ii] luminosity from z « 6 — 7 


LAEs and LBGs data are indicated with empty sym- 


bols ( 

Ouchi et al. 112013 Ota et al 

120141 Gonzalez-Lopez 

et al. 

2014 

Schaerer et al. 

2015 

Maiolmo et al. 2015). 

The recent |C ii| detection 

in the vicinity of a z « 7.1 


while 


LAE (Maiolino et al. 2015) is p lotted as a fil l ed sq uare, 
C III data at z ' 


6 by Capak et al. (2015) and 


Willott et al. (2015) are indicated with hlled hexagons 
and triangles, respectively. The [C ii]-SER relation pre¬ 
dicted by our model fairly reproduces the slope of the 
relation found in local dwarfs, as well as its trend with 
metallicity, although the scatter in the data is large. As 
shown in Eig. 5, the [CII] emission arising from the dif¬ 
fuse medium is expected to be < 40% in local galaxies, 
where the CMB attenuation on the CNM luminosity is 
negligible. This implies that the [CIIJ-SER relation is al¬ 
ways driven by the correlation between the SER and the 
intensity of the [CII] emission arising from PDRs, with 
and without taking into account CMB effects. Possible 
variations in the diffuse medium contribution to the total 
[CII] emission may result in a slight tilt of the [CII]-SFR 
relation, certainly within the current la scatter. 


4. SUMMARY AND DISCUSSION 

By coupling radiative transfer cosmological simulations 
of a z = 6.6 galaxy with a sub-grid ISM model and a PDR 
code (UCL_PDR), we have computed the [C ii] emission 
arising from the diffuse cold neutral medium (CNM) and 
molecular clouds in early galaxies, characterized by SFRs 
ranging from 0.1 to 100 MQyr“^. We have distributed 
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regions might have either a positive effect, by trigger- 
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Fig. 6.— [C Il]-SFR relation for the Cl, C02, and C005 models 
as a function of the slope (N) of the KS relation, with N ranging 
from 0.63 to 2.0. We indicate with thin dotted lines the [C Il]-SFR 
relation obtained for our fiducial model with {N = 1). The [C ll] 
luminosity arising from all gas phas es, and from molecular phase 
only, predicted byjOlsen et al.|l|2015|l for a sample of z Fa 2 galaxies 
with metallicities ranging trom U.4 Ziq to 1.67 Zq is plotted with 
green (purple) dot-dashed lines respectively. Green empty (filled) 
points represent upper limits (detections) of [C ll] in LAEs and 
LBGs at z Fa 5 — 7, as in Fig. 

metals in the ISM, both uniformly and accordin g to the 
Z — A relation found by Pallottini et al. (2014aI, to sim¬ 
ulate gas metallicities in the range U.U5-1 Zq. 

We find that the [C ii] line from high- 2 ; galaxies is dom¬ 
inated by emission from PDRs, while the CNM accounts 
for < 10% of the total flux. This is due to the fact that 
at these early epochs the CMB temperature approaches 
the spin temperature of the [C ii] transition in the CNM 
(Tcmb yCNM ^ 20 K) suppressing the flux contrast. 
The [C ii] spectrum predicted by our model is complex. 
It shows a pronounced peak (FWHM~ 50 km s“^) due 
to centrally located {v = 0) PDRs, and weaker [C ii] dis¬ 
placed {v ^ 200 km s“^) peaks from MCs in the galaxy 
outskirts. 

The predicted [C iij-SFR relation reproduces the corre¬ 
sponding relation found in local dwarfs remarkably well. 
Current upper limits from observations of z ^ 6 — 7 
LAEs and LBGs seem to indicate that these galaxies are 
characterized by a [C ii] luminosity fainter than expected 
from the local relation. Although this conclusion is still 
not definitive, it must be noted that the SFRs quoted 
for high-z galaxies are inferred from observations of the 
Lya emission line, and therefore must be considered as 
lower limits to the actual value. This implies that green 
arrows in Fig. [^should be moved towards higher SFR 
values, hence exacerbating the inconsistency with the 
local relation. Our results contain a caveat: we 

have so far neglected the possible effect of stellar feed¬ 
back (i.e. photo-evaporation, radiation pressure, [H ii ] 
thermal pressure) on molecular clouds. Broadly speak¬ 
ing, these effe cts should act to reduce the mass of the 
molecular g as (Tasker & Tan||2009 Tasker 2011 Tasker 
et al.|[20l5). On the other hand, the expansion of [H ii J 


et al. 

2011 

Haworth & Harrie 

2012 ), or disperse the 

surrounding cloud (Dale et al. 

2005). 'I'ypicai negative 


et al. 2006 Walch et al. 


2012 ), namely the age of stars 


taken into account m our radiative transfer calculations. 
This would imply a steeper scaling between the SFR and 
molecular hydrogen surface densities, because for a given 
value of SFR the mass of H 2 is lower. 

In Fig. we plot the [C ii]-SFR relation for models 
with constant metallicity (Cl, C02, and C005) as a func- 
tion of the slope in the range 0.63 < < 2.0 (iShetty 

| 2012 |. 




et al. 2013 


Bigiel et al(l[2008 
We show for reference the [C ii 
for our fiducial model A^ = 1 with thin dotted lines. 


Narayanan et 
Sh'R relation obtained 

At 


a given Z the steeper is the slope of the KS relation, the 
shallower is the [C ii]-SFR curve. Hence, a [C ii] deficit 
in z ^ 6 — 7 galaxies, if confirmed by deeper observa¬ 
tions, would favor a scenario in which star formation in 
early galaxies blows the molecular gas apart, reducing 
the amount of material from which most of the [C ii] 
emission arises. The deviation from the local [C ii]-SFR 
would then imply a modified Kennicutt-Schmidt relation 
in z > 6 galaxies. Stellar feedback effects are likely to be 
stronger in regions of very active star formation, more 
often located in galactic centers. This is particularly im¬ 
portant in high-z galaxies that are known to be more 
compact than their low-z counterparts. If so, negative 
feedback should preferentially suppress the peak in the 
[C ii] spectrum at the systemic redshift of the galax. 
Alternatively/in addition, as can be noted in Fig. 
the deficit might be explained by lower gas metallicities. 
In the same Figure, we also test our results with the 
[C ii]-SFR relation for all gas phases (PDR -|- molecu¬ 
lar -|- ionized) and from molecular phase only, calculated 
by Olsen et al. (2015) in a sample of simulated z « 2 
galaxies. As they point out in their work, their PDR 
component is best identified with what we call CNM, 
while their molecular component, located at the center, 
is comparable to what we call PDR emission. We find 
a nice agreement between our results and their findings 
for the molecular gas. This is something expected, given 
that the fraction of [C ii] emission arising from the CNM 
is almost totally attenuated due to the increased CMB 
temperature at z 6 . 6 . 

Finally, we note that the MC density distribution 
may play a role. Our simulated galaxy is character¬ 
ized by a mean molecular hydrogen number density 
rici ~ lO^’^cm”^; this quantity depends on the square of 
the assumed Mach number A4 = 10. Calculations per¬ 
formed with UCL_PDR show that molecular clouds char¬ 
acterized by densities 10 times higher (lower), for a fixed 
gas metallicity (e.g. log(Z/Z 0 ) = —1.5), would result 
into a [C ii] emissivity 5 times higher (20 times lower) 
than found here. Although we consider such large varia¬ 
tions of the Mach number unlikely, at present we cannot 
exclude that the corresponding shift in the mean MC 
density plays some role in the interpretation of the re¬ 
sults. 
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TABLE 1 

Summary of the models considered and relative results. 


Name 

( — ) 

profile 


FdiffIFtot (%) 


Lcil (10® Lq) 


1 M 0 yr 1 

10 Moyr-l 

100 M© yr-l 

1 M© yr 1 

10 M©yr-i 

100 M© yr-l 

Cl 

1 

cnst 

44(8) 

37(10) 

4(0.6) 

0.1 

1.6 

21 

C02 

0.2 

cnst 

26(4) 

20(5) 

2 ( 0 . 2 ) 

0.05 

0.7 

9.0 

C005 

0.05 

cnst 

41(8) 

20(5) 

3(0.2) 

0.007 

0.09 

1 

P005 

0.05 

A-Z 

24(4) 

21(5) 

2(0.3) 

0.05 

0.6 

8.0 


Note. —Name: model name, (Z/Zq): mean metallicity in solar units, profile: type of metallicity profile adopted, percentage 

of the [C ll] emission arising from the CNM without (and with, in bold) the attenuation due to the increased CMB temperature, Lqh: 
predicted [C ll] luminosity in 10* Lq. 
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